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Abstract 

The theoretical formulation of x-ray resonant magnetic scattering from rough surfaces and in- 
terfaces is given for specular reflectivity. A general expression is derived for both structurally and 
magnetically rough interfaces in the distorted-wave Born approximation (DWBA) as the frame- 
work of the theory. For this purpose, we have defined a "structural" and a "magnetic" interface 
to represent the actual interfaces. A generalization of the well-known Nevot-Croce formula for 
specular reflectivity is obtained for the case of a single rough magnetic interface using the self- 
consistent method. Finally, the results are generalized to the case of multiple interfaces, as in the 
case of thin films or multilayers. Theoretical calculations for each of the cases are illustrated with 
numerical examples and compared with experimental results of magnetic reflectivity from a Gd/Fe 
multilayer. 



I. INTRODUCTION 



X-ray reflectivity and offspecular diffuse scattering methods have been widely apphed 
over the last decade to characterize the morphology of rough surfaces and interfaces, partic- 
ularly with the availability of sources of ever-increasing brilliance for x-ray radiation. Similar 
techniques using neutron beams have also become widespread, particularly for the study of 
magnetic multilayers. In the case of x-rays, however, clement-specific information regarding 
the magnetic structure can be readily obtained by tuning the photon energy to that of an 
L-edge (in the case of transition or rare-earth metals) ^'^ or of an M-edge (in the case of 
actinides).^'^ The resonant enhancement of the scattering by magnetic atoms at such ener- 
gies can result in a large enough signal to be comparable to the dominant charge scattering. 
Resonant x-ray scattering at the K-edges of transition metals^ has also been used to obtain 
information about the magnetic structure, although the enhancement is not as large. Reso- 
nant magnetic scattering corresponds to the real part of the scattering amplitude, while the 
(absorptive) imaginary part gives rise to x-ray magnetic circular dichroism (XMCD), which 
has been used to obtain the values of spin and orbital moments in ferromagnetic materi- 
als. Detailed descriptions of the formalism for the interaction of x-rays with magnetically 
polarized atoms have been given in the literature,^"^'^ from which a complete description of 
magneto-optic phenomena in the x-ray region can be obtained and applied. 

Several resonant x-ray specular reflectivity experiments have been performed to obtain 
the magnetization within the layers of magnetic multilayers.^ '^^"^^ The analysis of these re- 
sults has generally used recursive matrix techniques developed for magneto-optics in the 
case of resonant x-ray reflectivity.^^ In general, roughness at the interfaces has been ignored 
or taken into account in an ad-hoc manner. In principle, representing roughness in terms 
of a graded magnetization at the interface and using slicing methods could enable one to 
calculate the effect of magnetic roughness on specular reflectivity at the expense of con- 
siderable computational effort. Rohlsberger has developed a matrix formalism (originally 
developed for nuclear resonant x-ray reflectivity) from which specular reflectivity incorpo- 
rating roughness can be calculated. It was not considered in his paper, however, that the 
magnetic interfaces can have different roughnesses from the structural (chemical) ones. In 
this paper, we define separately a structural and a magnetic interface to represent the actual 
interfaces and present analytical formulae taking into account both interface roughnesses. 
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which provide much faster computational method than the shcing methods and show good 
agreement with estabhshed formulae for chemical interface roughness. 

Methods were developed earher to calculate analytically the specular component of the 
charge scattering of x-rays by rough surfaces and interfaces using the Born approximation 
(BA) and the distorted- wave Born approximation (DWBA).^'''^^ The BA results were ex- 
tended to magnetic interfaces in an earlier publication^^ and have already been applied 
to interpreting x-ray resonant magnetic specular reflectivity measurements from magnetic 
multilayers.^^ However, the BA or the kinematical approximation breaks down in the vicin- 
ity of the critical angle and below, since it neglects the x-ray refraction. On the other hand, 
the DWBA takes account of dynamical effects, such as multiple scattering and the x-ray re- 
fraction, which become significant for smaller angles close to the critical angle and even for 
greater angles at the resonant energies or with soft x-rays. We present here the generaliza- 
tion of the DWBA to the case of resonant magnetic x-ray reflectivity from rough magnetic 
surfaces or interfaces. The principal complication is, however, that we now have to deal 
with a tensor (rather than scalar) scattering length, or equivalently an anisotropic refractive 
index for x-rays. This leads in general to two transmitted and two reflected waves at each 
interface for arbitrary polarization, which complicates the DWBA formalism. 

The plan of this paper is as follows. In Sec. II, we discuss a simple conceptual model 
for a magnetic interface and its relationship to the chemical (i.e., structural) interface and 
define the appropriate magnetic roughness parameters. In Sec. Ill, we discuss the (known) 
scattering amplitudes for resonant x-ray scattering and their relationship to the dielectric 
susceptibility to be used in the DWBA. In Sec. IV, we present the derivation of the scattering 
in the DWBA for a single interface with both structural and magnetic roughnesses. In Sees. 
V and VI, we derive the formulae for specular reflectivity from a magnetic interface using 
the self-consistent method in the framework of the DWBA and discuss numerical results. 
Finally, in Sees. VII- IX, we discuss the extension of the formahsm to the case of the specular 
reflectivity from magnetic multilayers and present some numerical results with experimental 
data from a Gd/Fe multilayer. In the following paper, we derive the formulae for the 
diffuse (off-specular) scattering from magnetic interfaces in both the BA and the DWBA. 
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II. MODEL FOR MAGNETIC INTERFACE 



Consider an interface between a ferromagnetic medium and a nonmagnetic medium 
(which could also be free space). Due to the roughness of this interface, the magnetic 
moments near the interface will find themselves in anisotropy and exhange fields, which 
fluctuate spatially (see Fig. 1). 

This will produce disorder relative to the preferred ferromagnetic alignment within the 
magnetic medium. A similar situation can arise at an interface between a ferromagnetic 
medium (FM) and an antiferromagnetic medium (AFM), where there is a strong antiferro- 
magnetic coupling between spins in the FM and the AFM. Random steps will then produce 
frustration in the vicinity of the interface, resulting in random disordering of the magnetic 
moments near the interface. Clearly in general correlation will exist between the height 
fluctuations of the chemical interface and the fluctuations of the spins, but a quantitative 
formalism to account for this in detail has not yet been developed. We make here the sim- 
plifying assumption that the ferromagnetic moments near the interface (or at least their 
components in the direction of the ferromagnetic moments deep within the FM layer, i.e., 
the direction of average magnetization M) are cut off at a mathematical interface, which we 
call the magnetic interface and which may not coincide with the chemical interface, either 
in its height fluctuations or over its average position, e.g., if a magnetic "dead layer" exists 
between the two interfaces (see Fig. 1). The disorder near the interface is thus represented 
by height fluctuations of this magnetic interface. The basis for this assumption, which is 
admittedly crude, is that the short (i.e., atomic) length-scale fluctuations of the moments 
away from the direction of the average magnetization give rise to diffuse scattering at fairly 
large scattering wave vectors, whereas we are dealing here with scattering at a small wave 
vector q, which represent the relatively slow variations of the average magnetization density. 
The actual interface can be then considered as really composed of two interfaces, a chemi- 
cal interface and a magnetic interface, each with their own average height, roughness, and 
correlation length, and, importantly, in general possessing correlated height fluctuations. 
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III. RESONANT MAGNETIC X-RAY SCATTERING AMPLITUDE 



The amplitude for resonant magnetic scattering of x-rays has been derived by Hannon 
et al.,^ and a discussion of the general formalism may be found in the review by Hill and 
McMorrow.^ There are two cases of practical importance, namely dipole and quadrupole 
resonances. We shall restrict ourselves here to the most commonly used dipole resonance, 
which is related to the L-cdges of transition metals and rare-earth atoms. The tensor 
amplitude for scattering /^^ from a magnetic atom is given by 



OTT 

.3A 



(e} • e,) 



- ^— (Fn-Fi_i)(e}xe,)-M 
+ ^(2Fio - Fn - i^i-i)(e} • m) (e, ■ m), (3.1) 

where e^, e/ are, respectively, the unit photon polarization vectors for the incident and 
scattered waves, M is a unit vector in the direction of the magnetic moment of the atom, 
A is the x-ray photon wavelength, /o is the usual Thomson (charge) scattering amplitude 
[/o = —ro{Z + f — if")], where tq is the Thomson scattering length (e^/mc^), Z is the 
atomic number, /'(< 0) and /"(> 0) are the real and imaginary non-resonant dispersion 
corrections. Flm is the resonant scattering amplitude, as defined in Ref. 6, and has the 
resonant denominator Ej.es — E — ir/2, which provides the resonance when the photon energy 
E is tuned to the resonant energy E^es close to the absorption edges. The lifetime of the 
resonance F is typically 1 — 10 eV, so that the necessary energy resolution is easily achivable 
at synchrotron radiation beamlines. (We assumed that q, the wave-vector transfer, is small 
enough here that the atomic form factor can be taken as unity.) Equation (3.1) has both 
real and imaginary (i.e., absorptive) components. The latter gives rise to the well-known 
phenomenon of x-ray magnetic circular or linear dichroism, whereas the real part gives rise 
to the scattering. Equation (3.1) yields 

/a/3 = AS^p -iBY, e^p^M^ + CM^M^, (3.2) 

7 

where 

A = /o + ^(i^ii + i^i-i), 
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C = — (2Fio-Fn-Fi_i), (3.3) 

OTT 

and a, /9 denote Cartesian components, and CaP'y is the antisymmetric Levi-Civita symbol 

{^xyz ~ ^yzx ~ ^zxy ~ 1) ^xzy ~ ^yxz ~ ^zyx ~ 1; Othcr C^p^ — 0). The diclcctric 

susceptibility of a resonant magnetic medium is given by 

^resonant(r) ^ ^^(r)/,^(r), (3.4) 

where /cq — 27r/A, n^(r) is the local number density of resonant magnetic atoms, and the 
variation of fap{r) with r reflects the possible positional dependence of the direction of 
magnetization M. The total dielectric susceptibihty is given by 

47r r f ~i 

Kq L 

- iBumiv) J2 eapjM^ir) + Cn^(r)M«(r)M^(r)j , (3.5) 

where po(r) represents the electron number density arising from all the other nonresonant 
atoms in the medium modified by their anomalous dispersion corrections when necessary. 
Using the constitutive relationship between the local dielectric constant tensor ea/3{r) and 
Xa/3(r), 

ea/3(r) = Sai3 + Xa/sir) ■ (3.6) 

We note that the magnetization gives the dielectric tensor the same symmetry as in conven- 
tional magneto-optic theory, namely an antisymmetric component linear in the magnetiza- 
tion. 



IV. THE DISTORTED- WAVE BORN APPROXIMATION FOR A SINGLE MAG- 
NETIC INTERFACE 

The results for specular reflectivity in the Born approximation (BA) have been derived 
in Ref. 19 and will be also summarized briefly in connection with the cross section in 
the following paper. ^'^ Here we discuss the scattering in terms of the distorted-wave Born 
approximation (DWBA). While this is more comphcated algebraically, it provides a better 
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description than the simple kinematical approximation or BA in the vicinity of regions where 
total reflection or Bragg scattering occurs. This treatment is a generalization of that used 
in Ref. 17 for charge scattering. The wave equation for electromagnetic waves propagating 
in an anisotropic medium with a dielectric susceptibility tensor given by Eq. (3.5) may be 
written as 

J2 [( V + kl)6ap - VaVp + klxap] Ep{r) = 0, {a, (3 = x,y,z), (4.1) 
P 

where E(r) is the electric field vector. 

Consider a wave incident, as in Fig. 1 with wave vector kj in the {x, z) plane {ki^y = 0) and 
polarization (/x = cr or tt), from a nonmagnetic (isotropic) medium for which Xap — Xo^o-P 
onto a smooth interface at z = with a magnetic medium, for which Xap is constant for 
^ < 0. 

Let us write for 2; < 

where the term x^ap P^^^ ^^^^ specifically depends on the magnetization M, as defined 

in Eq. (3.5). The incident wave (chosen for convenience with unit amplitude) may be written 
as 

= e.e^"^-'-. (4.3) 

This incident wave will in general give rise to two specularly reflected waves (where the 
index // refers to o" or tt polarization) and two transmitted (refracted) waves in the magnetic 
medium. The complete solution for the electric field in the case of the smooth magnetic 
interface is then given by 

E(k„^)(r) = e^e^^--+ E i?S(k,)e,e^^r-, z > 0, 



= E ^<0, (4.4) 

where k[ is the specularly reflected wave vector in the nonmagnetic medium, i/ denotes 
the polarization of the appropriate reflected component, the index j{— 1, 2) deflnes the 
component of the transmitted wave in the magnetic resonant medium with polarization 
{^j=i,2 = e^-*^^ and e^^^ respectively, as deflned in Appendix A), and k*(_7) the appropriate 



wave vector for that transmitted wave. The polarization vectors e may be real or complex 
allowing for linear or elliptically polarized waves. We denote such states in Eq. (4.4) 
quantum-mechanically by |kj,/x >. 

R^^^ and Tj'^ denote the appropriate reflection and transmission coefficients for the 
smooth surface and are expressed in terms of 2 x 2 matrices using the polarization bases for 
the incident and reflected (or transmitted) waves. The polarization basis is given by {e^, 
e^r), as shown in Fig. 1, for the waves in the nonmagnetic medium and {e^^\ e*^^^), as defined 
in Appendix A, for those in the magnetic resonant medium, respectively. The convention 
in which the polarization state of the refiected (or transmitted) wave precedes that of the 
incident wave is used for the subscripts in R^^} and T^^} , and the Greek and Roman letters 
are used for the polarization states in the nonmagnetic and magnetic medium, respectively. 
The explicit expressions of R!>^^ and T^^ for small angles of incidence and small amplitudes of 
the dielectric susceptibility and for special directions of the polarization and magnetization 
(i.e., M II X as shown in Fig. 1) are given in Appendix A. 

We should mention, however, that these specific conditions considered in Appendix A 
(and also in all other appendices) are reasonably satisfied for hard- and medium-energy x- 
rays and also for soft x-rays around transition-metal L-cdgcs with small angles (i.e., when 
01 -C 1 for the incidence angle Qj). We should also mention that, even when M is not 
parallel to the x-axis in Fig. 1, the expressions derived in the appendices can be still 
applied by considering only the x-component of the magnetization vector M. This is because 
the y- and z-components of M contribute negligibly to the scattering in comparison with 
with the dominant factor B = (3A/87r)(Fii — Fi_i) in Eq. (3.2) at small angles^^ when 
— ^ |2Fio — Fxi — -Fi-ij, which is generally satisfied for transition-metal and 

rare-earth L-edges.^ 

We note that the continuity of the fields parallel to the interface requires that 



where ()|| denotes the vector component parallel to the interface. 

We now discuss the structurally and magnetically rough interface. For this purpose we 
shall assume that the average height (along z) of the structural and magnetic interfaces is 
the same, i.e., we ignore the presence of a magnetic dead layer. This may be treated within 
the DWBA as simply another nonmagnetic layer and thus discussed within the formalism 




(4.5) 
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for treating multilayers as discussed in Section VII. We can write 

3(r) + A^^(r) + A-^ 



Xap{r) = xS(r) + AUr) + A-,(r), (4.6) 



where 



:2) 



= XlSa^ + A z<0, (4.7) 



^a/3(r) = (xi-Xo)5a/3, for < z < 5zc{x,y) if 5^c(a;,|/)>0 
= -(Xi -Xo)5a/3, for Szc{x,y) < z <Q if 52;c(a;,y)<0 
= elsewhere, (4.8) 



and 



= Xal for 0< z<5zmix,y) if Szmix,y)>0 

= -xi% for 52;m(a;,y) < ^ < if 6z,n{x,y)<0 

— elsewhere, (4.9) 

Szc{x, y) and 52;m(a^, u) define the structural (chemical) and magnetic interfaces, respectively. 

We may also define the time-reversed function corresponding to a wave incident on the 
interface with vector (— ky) and polarization u as 



A=<T,7r 

= E *(-k/)e,e^'^*/(^>, z<0, (4.10) 

where (— k^) is the wave vector of the wave specularly reflected from (— kj), and kj(j)^ 
is the wave vector of one of the two transmitted waves in the medium emanating from 
(— kj) incident on the surface, as shown in Fig. 2. Note that, for consistency with the 
conventions used in Eq. (4.4), the polarization vectors in Eq. (4.10) are defined in the 
ordinary coordinate system where their phases are considered along the left-to- right direction 
in Fig. 1. Otherwise, the polarization vectors in Eq. (4.10) should be replaced by their 
complex conjugates. 

We have also the conditions 

(k,)|| = (k;^)ii = (k*(j)),,. (4.11) 



The DWBA then yields the differential cross section for scattering by the rough interface 
from (kj, n) to (k/, u) as 

^^^(iTffV (4.12) 



, /il2 

dn 167r2\'' 

where T-^^ —< k/, i/|T|kj, > is the scattering matrix element, and ^...^ in Eq. (4.12) 
denotes a statistical averaging over random fluctuations at the interface. Following Ref. 17, 
we split the cross section into two parts: 



da 

dn 



167r2 



+ 



167r2 



(4.13) 



The first term in Eq. (4.13) represents the coherent (specular) part of the scattering, which 
corresponds to a statistical averaging of the scattering amplitude, and the second term 
corresponds to the incoherent (diffuse) scattering. In this paper, we shall deal with the first 
term only, while the diffuse scattering will be addressed in the following paper. ^° 

The DWBA consists of approximating the scattering matrix element by the expression 

< k^,i.|r|k,,/. > = ki < -^Mx^'^Ki^) > 

+ kl< -kj, u\A'\ki, n>^kl< -kj, i/| A™|ki, ji > . (4.14) 

Here |E^(r) > denotes the "pure" incoming wave in Eq. (4.3), | — kj, v > denotes the state 
in Eq. (4.10), and the matrix clement involves dot products of the tensor operators X^^\ ^'^^ 
and A™ with the vector fields < — kj, v] and < kj, While x^^^ represents an ideal system 
with a smooth interface, A'^ and A"* are perturbations on due to interface roughnesses. 

For the smooth surface, only the first tensor is nonvanishing, and, following Ref. 17, we 
can show from Eqs. (4.3) and (4.10) that 



kl < -kJ,i/|x^°^|E;(r) > = iAkl5k,^kfJk,ykfy 



j a/3 

X r dze-'^^^f^^^-^^'^^' , 

J —oo 



(4.15) 



where A is the illuminated surface area, and -R^'^^(kj) is the reflection coefficient for the 
smooth surface, as defined in Eq. (4.4). The details of Eq. (4.15) are presented in Appendix 
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B. By comparison with Eq. (4.15) for the smooth surface, the scattering matrix element for 
the rough surface in Eq. (4.14) can be analogously defined by 

< k/, i^|T|ki, n >= 2iAki,R^^{ki)5ki^kfJkiykfy, (4.16) 

where i?i,^(kj) denotes the reflection coefficient for the rough surface. 

On the other hand, for the reverse case where a wave is incident from a resonant magnetic 
medium to a nonmagnetic (isotropic) medium, similarly to Eq. (4.15), the scattering matrix 
element for the smooth surface can be shown to be 

kl < -kj,/|x('^)|k„i >= AiAkUj)R^r]{h)Sk,.kJk,yk,y: (4.17) 

where the incoming wave from the resonant magnetic medium |kj,j > is used instead of 
the "pure" incoming wave from the vacuum E^r) in Eq. (4.3). The use of Eqs. (4.15) 
and (4.17) in Eqs. (4.14) and (4.12) in the case of the smooth surface and the derivation 
of the corresponding reflectivity in the usual manner, as discussed in Ref. 17, shows that 
Eqs. (4.15) and (4.17) must be identically true. Similarly to Eqs. (4.15) and (4.16), the 
scattering matrix element for the rough surface between reversed layers can be also defined 
by analogy from Eq. (4.17) as 

< k/, j'|T|ki, j >= 4iAki^{j)Rj,j{ki)6k^^kfJkiykfy, (4.18) 
where Rj/j(ki) denotes the reflection coefficient for the rough surface between reversed layers. 

V. REFLECTION AND TRANSMISSION COEFFICIENTS USING THE SELF- 
CONSISTENT METHOD 

To calculate specular reflectivity, we make an approximation in the spirit of Nevot and 
Croce.^^ To evaluate the matrix elements in Eq. (4.14) involving A^^ and A^^, we assume 
for E(kj,/i) in Eq. (4.4) the functional form for z > analytically continued for z < 0, 
while for the time-rcvcrscd state E"^(— kj,z/) in Eq. (4.10) the functional form for 2; < 
analytically continued to 2; > 0. Then, bearing in mind that for specular reflectivity ky = k[ 
and using Eq. (4.5), we obtain for the statistically averaged amplitude (^■^*): 

(kl < -kJ,i.|A'='-|k„/. >) = zAkl ?i?(-k/) 

11 



X 



* \c,m 



a/3 



g-i9uO)52c,m(a:,J/)\ _ ^ 



S-a^a/3 eA/3 



A=iT,7r 



a/3 



-iq2z{j)Szc,m{x,y) 



- 1 



(5.1) 



where 



(5.2) 



and A^'J is the value defined for < 2; < 5zc,m in Eqs. (4.8) and (4.9). Prom Eqs. (4.15)- 
(4.16) and (5.1), we see that, at the specular condition, we can write Eq. (4.14) as 



(5.3) 



where 



rr - sr kl 



(Xi-Xo)Ee-ae.a[e-^''-(^>— 1] 



al3 



(5.4) 



and replacing qiz, in by q2z, e.\ produces V^^a- Here we made the customary Gaussian 
approximation for the height fluctuations Szc,m{x,y), and ctc, am are the root-mean-squared 
structural and magnetic roughnesses, respectively. Note that the correlation term C/,^^ due 
to the roughness in the reflection coefficient contains only independent contributions of 
chemical and magnetic roughnesses expressed via Uc and a^, respectively. According to Eq. 
(4.13), the diffuse scattering must contain the cross-correlation component due to the term 

(ir/f). 

A better approximation than Eq. (5.3) may be obtained by using the rough-interface 
refiection coefficient R^n instead of the smooth-interface R^^ in the wave functions of Eqs. 
(4.4) and (4.10), thus getting a self-consistent matrix equation in terms of the 2x2 matrices, 
R, U, V. This leads to 



R = R(°) + U + VR, 



(5.5) 



whose solution is 



R 



V)-i(R(°) + U). 



(5.6) 
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Similarly, for the reverse interface between upper resonant magnetic and lower nonmagnetic 
layers, we can have the same solution as Eq. (5.6) from Eqs. (4.17) and (4.18). The exphcit 
expressions of U, V, R^'^^ matrices in Eq. (5.6) for both cases are given in Appendix C. 

For nonmagnetic interfaces, the matrices are all diagonal((T and tt polarizations are de- 
coupled), and it has been shown that Eq. (5.6) leads to the familiar Nevot-Croce form^^ for 
the reflection coefficient, i.e.. 

The derivation of this is shown in Appendix D. For the magnetic interface, this simplified 
form for the refiection coefficient does not have any analogue. Nevertheless, at sufficiently 
large values of the reflectivity takes the familiar Gaussian form i?(°)e~^^'^eff. However, a^Q 
does not always take the form predicted by the simple kinematical theory [i.e., cr^ for cr — > cr 
reflectivity, for cr — > tt reflectivity, and |(c^ + cr^) for (7_|. — 7_) in the case of circularly 
polarized x-rays] as we shall see in the numerical example shown below, which provides a 
counter-illustration of the rule that, at large q^, the DWBA becomes identical to the Born 
approximation or kinematical limit. 

For circularly polarized incident x-rays with e±(ki) = (ea(ki) ± ieT^{ki)^ /V^, the reflec- 
tion amplitudes for a- and 7r-polarization are given by 

("A^J^-]. (5.8) 

where R is the 2x2 matrix reflection coefficient in Eq. (5.6) . The reflected intensities without 

polarization analysis for the outgoing beam, / = \J\R^^^~+Jr^ , can be then evaluated for 
the opposite helicities of incident beams as 

1+ - I- = 2 lm[RnRl2 + R21RI2], (5-9) 

where Rij is the ij-element of the 2x2 matrix R. 

Since Parratt's recursive formula for multiple interfaces includes only reflection coeffi- 
cients, its extension to the rough interface case does not need the transmission coefficient 
to account for interface roughness. On the other hand, in our case where the fields are not 
scalars, the transmission coefficients are requisite to calculate recursive 2x2 matrix formu- 
lae for multiple magnetic interfaces, which will be discussed in Sec. VII. For completeness. 
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therefore, let us now calculate the transmission coefficient Tj^ from a rough interface. In the 
spirit of Ref. 22, we assume for E(kj,//) and E^(— k/,,;) the functional forms analytically 
continued both for 2; > and for z < as follows: 

E(k„/.)= J: Tj?lik,)ej,e'^mr^ (5.IO) 
E^(-k/,j)= E Tif (-k,)e.e-*/-, (5.11) 

where Tj^^*(— ky) in Eq. (5.11) denotes the transmission coefficient "from" a mag- 
netic (anisotropic) medium "to" a nonmagnetic (isotropic) one, whose explicit form is given in 
Appendix A. For the smooth surface, the scattering matrix element between the eigenstates 
I — kj, j > and |kj, n > can be then written as 



f j' a/3 

J— 00 

= AiAkUj)Tjl\k,)6,,^,J,,^,^^, (5.12) 

where 7^-°'*(kj) is the transmission coefficient for the smooth surface, as defined in Eq. (4.4). 
The details of Eqs. (5.12) are given in Appendix B. 

In comparison with Eq. (5.12) for the smooth surface, the scattering matrix element for 
the rough surface, as shown in Eq. (4.14), can be analogously defined by 

(5.13) 

where Tjn(ki) denotes the transmission coefficient for the rough surface. 
For the statistically averaged amplitude (T^^'), we obtain 



^-i<l3zU')^Zc,m{x,y)\ — 14) 



kl<-k^f,j\A'^'^\k,,^,>) = iAk',Y.TS\-kf) 

V 

if \c,m 
j' aP HZz\J ) 

and 

q3z{f)-^-kf.-kl{j'). (5.15) 
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Prom Eqs. (5.12)-(5.13) and (5.14), we see that we can write the scattering matrix 
element in the DWBA, as shown in Eq. (4.14), as 



/=1,2 



(5.16) 



where 



Tff(-k,) 



a 
a/3 



(5.17) 



In the same way as we did for the reflection coefficient, using the rough- interface trans- 
mission coefficient Tj^ instead of the smooth-interface Tj^ in the right side of Eq. (5.16), 
thus getting a self-consistent matrix equation in terms of the 2x2 matrices, T, V, gives 



whose solution is 



T = T(°) + V'T, 



(1_ V')-^T(°\ 



(5.18) 



(5.19) 



Similarly, for the reverse interface between upper resonant magnetic and lower nonmagnetic 
layers, we can also have the same solution as Eq. (5.19). The explicit expressions of V and 
T*^^'' matrices in Eq. (5.19) for both cases are given in Appendix C. 

For nonmagnetic interfaces, it is shown in Appendix D that Eq. (5.19) reduces to 



2^^2'(o)gi(|fc.|-|fe*|)'a?^ 
which has been found by Vidal and Vincent. 



(5.20) 



VI. NUMERICAL EXAMPLES FOR A SINGLE MAGNETIC SURFACE 

We now illustrate numerical examples of the above formulae calculated for a Gd surface 
with varying degrees of structural and magnetic roughness. We have considered only the 
case where the magnetization vector is aligned along the sample surface in the scattering 
plane in order to enhance the magnetic effect. 

Figure 3 shows the x-ray resonant magnetic reffectivities calculated at the Gd La-edge 
(7243 eV) from Gd surfaces with different interfacial widths for structural (uc) and magnetic 
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(am) interfaces. In Fig. 3(a)-(c), the interfacial width of the structural interface is larger 
than that of the magnetic interface, that is, ac — 8 A and (7^ = 3 A. On the other hand, 
in Fig. 3(d)-(f), the interfacial widths are reversed, that is, Uc = 3 A and — 8 A. In the 
kinematical approximation(BA) a — > a scattering (solid lines in the top panels of Fig. 3) 
corresponds to pure charge scattering, and a ^ ir scattering (dashed lines in the top panels 
of Fig. 3) to pure magnetic scattering, and the differences between the reflected intensities 
for right- (/+) and left- (/_) circularly polarized incident beams (circles in Fig. 3) correspond 
to the interferences between charge and magnetic scattering. 

Kinematically, the reflected intensities from each scattering channels are proportional 
to a simple Gaussian form, exp(— ci^g^), where a is the interfacial width of corresponding 
scattering channel, i.e., cTc for I^^a, <^m for Jo-^tt, and ^ [o1 + for (/_|_ — /_). The 

middle panel of Fig. 3 shows natural logarithms of the reflectivities from rough interfaces 
normalized to those from ideal systems without roughness as a function of the square of 
the wave vector, g^, whose slopes are then equal to the squares of the interfacial widths 
for their corresponding scattering channels. In Fig. 3(b), the slopes obtained from our 
dynamical calculation for the case of cxc = 8 A and cxm = 3 A show good agreement with 
the kinematical results mentioned above. On the other hand, in Fig. 3(e), the slopes of 
/ct-^tt a-nd (/+ — /_) for the opposite case, = 3 A and cr^ = 8 A, are not equal to the 
squares of their corresponding interfacial widths but follow the slope of /o— »ct a-t high q^^. 

This indicates that the kinematical argument mentioned above, i.e., one-to-one corre- 
spondence such as (7 — > TT channel to pure magnetic scattering, is no longer valid for such 
a case of larger magnetic interfacial width, as shown in Fig. 3(e). In other words, both 
contributions from charge and magnetic scattering should be taken into account for every 
scattering channel, which is naturally included in the dynamical theory (such as our self- 
consistent method). In the case shown in Fig. 3(e), since the charge-scattering channel is 
much stronger than the magnetic-scattering channel and also drops off much more slowly 
with qz due to decreased roughness, there is conversion of cr — > tt polarization at larger q^ 
even when the "pure" magnetic scattering has become negligible in the kinematical limit, 
because of magnetic scattering out of the still strong charge channel. Thus the o ^ -n and 
(/+ — /_) reflections will asymptotically decay at a rate governed by the decay of the charge 
channel, which is determined by Oc alone. 

However, it is not easy to flnd a physical system where a magnetic interfacial width is 
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larger than the structural one at the same interface, as shown in Fig. 3(f). Instead, such 
a rougher magnetic interface can occur in a magnetic system, where a magnetically "dead" 
layer exists near the top surface and so the average position of the magnetic interface may 
not coincide with that of the structural interface, as shown in Fig. 3(i). In Fig. 3(g) 
(J_i_ — /_) (circles) shows an oscillation due to a magnetically dead layer with its thickness 
of 20 A. In this case, the slopes in Fig. 3(h) follow again the kinematical result mentioned 
above because the magnetic interface and the structural one are separated spatially. 

As a further check on our calculations, we have calculated the reflectivity by dividing 
the error-function profile, as shown in the bottom panel of Fig. 3, into many very thin 
slices and using the 2x2 recursive matrix formulae without any roughness assumptions.^^ 
We found that the results using this slice method arc exactly the same as those from our self- 
consistent method assuming Gaussian height distributions in Fig. 3. Thus our self-consistent 
method based on the DWBA produces very accurate results for the x-ray resonant magnetic 
reflectivity and much faster computationally. 

VII. MULTIPLE MAGNETIC INTERFACES 

For a multilayer with multiple interfaces, each layer can be characterized by its dielectric 
susceptibility tensor Xai3,n for the n-th layer, which can be Xa/3,n = Xn^ap ioi nonmagnetic 
(isotropic) layers and Xai3,n — XnSai3+Xap,n fo^ magnetic (anisotropic) layers. For each rough 
interface, we can use the self-consistent DWBA to define the reflection and transmission 
coefficients, in the same way as in Sec. V, which are given by 

R„ = (I-V„)-i(Rf + U„) = K*, 

T„ = (I-Vg-iTW=M^*, (7.1) 

where R„, T„ are the reflection and transmission coefficients for the n-th rough interface, 
and R^°\ T^^ are those for the corresponding smooth interface. The explicit expressions 
for R^f), T^f), U„, V„, and V'„ matrices in Eq. (7.1) are given in Appendix C, depending 
on whether the upper and lower layers on the n-th interface are nonmagnetic or magnetic 
layers, respectively. 

By analogy with the recursion relation for the coupled waves derived for the smooth 
interfaces in Appendix E (originally developed by Stepanov and Sinha^^), introducing W^'^ 



17 



matrices for the rough interfaces, we may derive the recursion relation analogous to Eq. 
(E5), obtaining 

= ^„M-„ (7.2) 

where and are defined by 



n+l 

= (7.3) 

Finally, the specular reflectivity of a magnetic multilayer with rough interfaces can be 
obtained by 

i?o = W;jTo. (7.4) 

To calculate the sum and difference in the reflectivities for (+) and (— ) circularly polarized 
incident x-rays, substituting Tq — -^{1, ±i) in a similar way to Eqs. (5.8) and (5.9) yields 

/+ + /- = i(w^;,*)iip + i(w^;;)i2p + i(w^;,*)2ip + i(w^;,%r, 

= 2 im [{w;}Mw;,% + {w;}Mw;,%,] , (7.5) 

where {W]^)ij is the ij-element of the 2x2 matrix Wf^. 

The above suggested approach to calculating the effects of roughness in multilayers on 
specular reflectivity is an approximation analogous to those used previously in several publi- 
cations on charge-only roughness. Basically, it corresponds to averaging the reflection 
coefficient (or the scattering matrix) of each interface over the interface roughness. The 
comparison with the results of rigorous "slicing method" made in Ref. 26 has proven that 
such an approximation works very well. A possible reason for the excellent validity of this 
approximation is that the roughness effect is mainly displayed at greater incidence angles, 
where the reflection is small and the multiple scattering can be neglected (the total re- 
flection amplitude is a linear sum of contributions from individual interfaces). Note that, 
since we are considering the coherent scattering which involves only the statistical average 
of the scattering amplitude in Eq. (4.13), there is no contribution from any cross- interface 
correlations of roughness. This will not be the case with diffuse (off-specular) scattering. 
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VIII. NUMERICAL EXAMPLES FOR MULTIPLE INTERFACES 



We present here numerical examples for x-ray resonant magnetic reflectivity from a Gd/Fe 
multilayer using the above formulae. Since Gd/Fe multilayers (MLs) have vastly different 
Curie temperatures and strong interfacial coupling of Gd and Fe, these systems give rise 
to complex magnetic structures depending on the layer thickness, temperature, and applied 
magnetic field. Due to the advantage of Gd L— edge resonances available in the hard x-ray 
regime, several experimental studies from these Gd/Fe MLs have been performed using x-ray 
resonant magnetic reflectivity measurements. ^^'^^'^^ Again, we have considered only the case 
where the magnetization vector M || x. 

We have used the experimentally determined values for charge and magnetic resonant 
scattering amplitudes, /c,m. = /cm + ^/cm; ^^e resonant energy. The energy dependence 
of the absorption coefficient for opposite helicites, fi^{E), were measured from a [Gd(51 

A) /Fe(34 A)]i5 multilayer, which will be discussed below as an experimental example. The 
edge-step normalized /"^ were obtained from the charge and magnetic absorption coeffi- 
cients, iic,m W = (a*"^ + A«~)/2, iim^ IJ''^ - A«~], through the optical theorem, f"^ oc iic,m- 
Their absolute values were determined using the tabulated bare-atom scattering amplitudes 
away from resonance. Real parts were obtained from differential Kramers-Kronig transforms 
of imaginary parts. Figure 4(a) and (b) show the charge and magnetic scattering ampli- 
tudes around the Gd L2-edge obtained in such absorption measurements. These values are 
in good agreement with the calculated ones from the hsted values of A and B in Eq. (3.3) 
obtained from Ref. 8. For consistency of the definitions, it should be mentioned that the 

used here correspond to Im[A, B] in Eq. (3.3), whereas the /^^^ correspond to — Re[A, 

B] , respectively. 

Figure 5 shows the calculated x-ray resonant magnetic reflectivities from a [Gd(51 
A)/Fe(34 A)] 15 multilayer for different incident x-ray energies indicated in Fig. 4: (a) 
7926 eV, (b) 7929 eV, (c) 7931 eV, and (d) 7935 eV. The fines and symbols represent the 
sum and difference in the refiected intensities for (-I-) and (— ) circularly polarized incident 
x-rays, respectively, calculated using Eq. (7.5). Since the Gd/Fe multilayer was assumed to 
be sandwiched between Nb buffer (100 A) and cap (30 A) layers, the Kiessig fringes between 
the multilayer peaks in (/+ + /_) intensities result from the interference of the scattering 
of Nb layers and thus show little energy dependence around the Gd absorption edge. On 
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the other hand, (/+ — I-) intensities around the multilayer peaks show a clear energy de- 
pendence in signs and magnitudes relative to (/+ + /_) intensities. In Fig. 5(a) and (d) 
at which energies becomes much smaller than the signs and relative magnitudes of 
— /_) intensities follow simply the energy dependence of in Fig. 4(b), as expected 
in the kinematical approximation.^^ At the energies close to the absorption edge where 
cannot be neglected, however, one can hardly expect the signs and magnitudes of (/+ — /_) 
intensities to be obtained directly from the values of and in Fig. 4(b). Therefore, 
quantitative analysis on x-ray resonant magnetic reflectivity data at the resonant energy 
requires accurate calculation taking into account refraction and multiple scattering effects 
using dynamical theory, such as our self-consistent method presented above. 

In order to study the effect of the magnetic roughness amplitude, (/_|_ — J_) intensities 
for two cases, am < Oc and Um > cTc, have been calculated, as shown in Fig. 6. The 
calculations for — have been shown in Fig. 5. For all cases, the charge roughness 
amplitudes were assumed to be (7c,Fe/Gd = 4.7 Aand crc,Gd/Fe = 3.6 A. At the energy of 7935 
eV, the intensities of (/+ — /_) around the multilayer peaks are proportional to a simple 
Gaussian form, exp(— cr^g^), as shown in Figs. 5(d), 6(a), and 6(b). This is consistent 
with the kinematical calculations,^^ and a for (/_,. — /_) corresponds to \j (a^ + o"^)/2 as 
given by the kinematical argument. On the other hand, at the energy of 7929 eV where 

cannot be neglected, such a kinematical argument is no longer valid. Comparing Figs. 
5(b), 6(c), and 6(d), we can see that the magnitudes of (/+ — /_) peak intensities do not 
follow a Gaussian form, exp(— ct^q'^), but their signs change from negative (filled circles) 
to positive (open circles) values. This indicates that (/+ — /-), which is known to be the 
charge-magnetic interference scattering in the kinematical theory,^ is sensitive even to the 
interference between charge and magnetic roughness amplitudes. However, it should be 
mentioned again that this result cannot be reproduced by the kinematical calculation but 
only by the dynamical one presented above. 

Let us now consider the case where the magnetic structure in the resonant layers may 
not coincide with the chemical structure. For example, the ferromagnetic moments in Gd 
layers near Gd/Fe interfaces can be induced by the adjacent ferromagnetic Fe layers above 
the Curie temperature of Gd atoms, ^^'^^ or a magnetically "dead layer" may exist at an 
interface between a ferromagnetic layer and an antiferromagnetic layer. Here we assume 
simply three different magnetization depth profiles in the Gd layers of a Gd/Fe multilayer. 
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as shown in Fig. 7: uniform magnetization (A), ferromagnetic moments only near the Gd/Fe 
interfaces (B), ferromagnetic moments near the centers of Gd layers between magnetically 
dead layers (C). 

Figure 8 shows the results of calculations of x-ray resonant magnetic reflectivities from 
[Gd(51 A) /Fe(34 A)]i5 MLs with the different magnetic structures of Fig. 7. We assumed all 
magnetic roughness amplitudes of am — 4.2 A (effectively same as ctc) and the photon energy 
of £■ = 7929 eV. In Figs. 8(a)-(c), Gd layers were assumed to be magnetized only near the 
Gd/Fe interfaces [model (B)], and the thickness of each magnetized layer was assumed to 
be 4.6 A (a), 8.4 A (b), and 12.8 A (c). On the other hand, in Fig. 8(d)-(f), Gd layers 
were assumed to be magnetized in the middle of each Gd layer and sandwiched between 
magnetically dead layers [model (C)], and the thickness of each dead layer was assumed to 
be 4.6 A (d), 8.4 A (e), and 12.8 A (f). 

Unhke the case of uniform magnetization [model (A) in Fig. 7] shown in Fig. 5(b), 
(/+ — /_) intensities in Fig. 8 for models (B) and (C) show no suppression in peak intensities 
due to the charge-magnetic interference, as discussed above. This may be ascribed to a 
spatial separation between the charge and magnetic interfaces in models (B) and (C), as 
shown in Fig. 7. 

In addition, the signs and relative magnitudes of (/+ — /_) intensities at the multilayer 
peaks change remarkably as the thicknesses of magnetized layers change. In general, the 
peak intensities of the (m + n)— th order ML peak and its multiple orders are weak compared 
to other peak intensities when the thickness ratio between two constituent layers is n/m. 
For example, in our Gd(51 A)/Fe(34 A) multilayer, the fifth peak corresponds to such 
a suppressed peak. Therefore, different thicknesses of magnetic layers readily change the 
order of the suppressed peak in (/_|_ — /_) intensities, as shown in Fig. 8. On the other hand, 
the signs of (/+ — /_) intensities for models (B) [Fig. 8(a)-(c)] and (C) [Fig. 8(d)-(f)] are 
opposite each other, because their magnetic structures are exactly reversed. 

IX. EXPERIMENTS 

X-ray resonant magnetic reflectivities were measured from an Fe(34 A)/[Gd(51 A)/Fe(34 
A)] 15 multilayer. The multilayer was sputtered onto a Si substrate using Nb buffer (100 
A) and cap (30 A) layers. SQUID magnetometry and XMCD measurements show that the 
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multilayer couples antiferromagnetically at the Gd/Fe interfaces and have coercive fields 
< 50 Oe at 300 K. X-ray measurements were performed at sector 4 of the Advanced Photon 
Source at Argonne National Laboratory. Undulator radiation was monochromatized with 
double Si(lll) crystals and its polarization converted from linear to circular with a diamond 
(111) quarter-wave plate operated in Bragg transmission geometry. The sample was placed 
in a, B — 2.1 kG field parallel to its surface and in the scattering plane. Specular magnetic 
refiectivity was measured at room temperature with a photon energy near the Gd L2 res- 
onance (7929 eV) across multilayer Bragg peaks by switching the hehcity of the incident 
radiation at each scattering vector = (47r/A)sin^, with 9 being the grazing incidence 
angle. 

Figure 9 shows specular reflectivity curves obtained by adding [(a), (/+ -|- /_)] and 
subtracting [(b), (/+ — /_)] reflected intensities for opposite helicites of the incoming x- 
rays. Symbols represent measurements and solid lines represent the fits calculated us- 
ing Eq. (7.5). From the fit for (7+ -|- 7_) intensities, we obtained the layer thicknesses 
dcd — 50.74 ± 0.09 A and rfpe = 33.98 ± 0.09 A, and the roughness amplitudes of charge 
interfaces crc,Fc/Gd = 4.7 ± 0.1 A and (Jc.Gd/Fc = 3.6 ±0.1 A. From the fit for (/+ — /_) in- 
tensities, we found that the Gd layers were fully magnetized only near the Gd/Fe interfaces 
at room temperature, which is above the bulk Tc of Gd. This magnetization is induced by 
a strong antiferromagnetic exchange interaction with the magnetically ordered Fe layers. 
From the best fit, the thickness of the ferromagnetic Gd layer was estimated to be 4.5 ± 0.3 
A, which is consistent with our previous work.^^ Magnetic roughness amplitudes for Gd/Fe 
(Fe/Gd) and Gd-ferromagnetic/Gd-paramagnetic interfaces were estimated to be 4.2 ± 0.1 
A and 4.6 ± 0.1 A, respectively. 

X. CONCLUSIONS 

The formulae for x-ray resonant magnetic specular refiectivity have been derived for 
both single and multiple interfaces using the self-consistent method in the framework of the 
distorted- wave Born approximation (DWBA). For this purpose, we have defined a structural 

and a magnetic interface to represent the actual interfaces. The well-known Nevot-Croce 
expression for the x-ray specular reflectivity from a rough surface has been generalized and 
examined for the case of a magnetically rough surface. The formalism has been generalized 
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to the case of multiple interfaces, as in the case of thin films or multilayers. Numerical 
illustrations have been given for typical examples of each of these systems and compared 
with the experimental data from a Gd/Fe multilayer. We have also presented the exphcit 
expressions in the small-angle approximation, which are readily applicable to transition- 
metal and rare-earth L-edge resonant magnetic reflectivities. The code for the calculations 
in this paper is also available in C language by emailing to D.R.L. (drlee@aps.anl.gov). 
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APPENDIX A: EXPLICIT EXPRESSIONS FOR i^ffi, rj°^ USING 2x2 MATRIX 
FORMULAE 



To calculate the exphcit expressions for R^^l and Tj'^ in Eq. (4.4), we follow Stepanov and 



Sinha's approach^^ developed for magnetic resonant reflections from ideally smooth inter- 
faces. The electric field E2;<o(r) inside the magnetic medium with a dielectric susceptibihty 
tensor given by Eq. (4.2) can be represented as 

E^<o(r) = Ee-^*^°"^+''=o^°'^*'^, (Al) 

where 9i is the incidence angle, as shown in Fig. 1. The parameter u can be a complex 
number due to absorption or total reflection. Substituting this in the wave equation Eq. 
(4.1), we obtain 

J2 [(sin^ - u'^)Sap + naup + Xap] Ep = 0, (A2) 

where Ua — ka/ko, i.e., = cosOi, Uy — 0, and Uz — —u. 

If we consider the case where the magnetization vector is aligned along the sample surface 
in the scattering plane, i.e., M || x in Fig. 1, the tensor Xai3 of a resonant magnetic medium 
can be written from Eq. (3.5) as 

( Xi + C' Q \ 



'Mllx 



V 








Xi -iB' 
iB' Xi J 



, (A3) 
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where 



Xi = —r2Poir}ro + -r^Anm{r), B' = -^Bnm{f)M^, C = -^Cn^{f)Ml. (A4) 

Assuming that the incidence angle 6i is small (sin^, <S 1 and ria; = cos^j ~ 1) and even 
at the resonance Xap remain small (|Xa/3| ^1), and inserting Eq. (A3) into Eq. (A2), the 
dispersion equation for a nontrivial solution of Eq. (A2) can be then approximated by 



1 -u 

ef + xi-u^ -iB' 

-u iB' e? ■ 



Xi 



0, 



(A5) 



and the respective roots are ii(^'2.3,4) _ ±-^^1 + xi ± B'. Two roots of these w^-'^'s with 
Im[M(^'^'] > and the other two roots with Im[M*^^'^)] < correspond to transmitted and 
reflected waves in the medium, respectively. For each of the waves Eqs. (A2) and (A5) give 

(i = i,-,4) 



iB' 



iB' 



V ' 2/ o- ■ 



(A6) 



And if we denote 



u 



(1) 



ei + xi + B' = u+, u^^) = ^ef + xi-B' = u_, 



(A7) 



we may then write 

Ei'^ = iE^J\ eP = -zEj^\ Ef^ = zE(^\ ^(^) = -^Ei'\ 

EU) = „0)^0) (j = l,...,4). (A8) 

Since \u'^^'^ < 1, E^^'^ can be neglected, then the polarizations of the waves e^-') in the 
magnetic resonant medium can be reduced to the circular polarizations 



^^^^ ^ E^^^e^ + Ei^^e^ (y = e„ z^e^), 



(A9) 



If the wave field E2>o(r) with the incident and specularly refiected waves inside the non- 
magnetic (isotropic) medium can be represented as 



E,>o(r) = (Eoe-^'^"""^ + E^e^'^"""^) e"^ (uo = ^Jef + Xo) , 



(AlO) 
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the boundary conditions for the waves, E2>o(r) and E2<o(r) in Eqs. (Al) and (AlO) must 
be satisfied for the lateral components Ey and H|| of electric fields and magnetic fields, 
respectively. Since H oc [k x E] , this gives 



j 
j 

j 

Eo. + Er^ = J2(u^^^Ei^^ + n,Ei^^)^J2Ei'\ 

3 j 

where the approximation in the last equation was obtained by lu^-'^l <^ 1 and Ux 
Eqs. (A6)-(A8), the above equations can be expressed in the 4x4 matrix form 



(All) 



1. Using 



10 1 
10 1 

uq —uo 

Mo —Uq j 



\ f E ^ 



Eqtt 
ERa 
Ert: 



I 



1 



Ei^) 



(A12) 



111 

i —i i —i 

U+ U- — —U- 

y iuj^ —iu- —iu+ iu- j 

Representing the waves as the vectors Tq — {Eq^, Eq^^), Rq — {Er^, Er^^), Ti — {E^^\eP), 
and Ri = {E^\eI^^), the 4x4 matrices in Eq. (A12) can be reduced into four 2x2 blocks 



-to 
yRoj 



J 



Ti 



(A13) 



where X**, X**", X*"*, X^'^ can be obtained by multiplying the inverse of the 4x4 matrix at 
the left side of Eq. (A12) onto the both sides. Since the reflected waves inside the medium 
vanish for a single surface, E^^^ = E*^^) = [i.e., Ri — (0, 0)], the "unknown" waves Rq and 
Ti in Eq. (A13) can be expressed via the "known" waves Tq and Ri as 




f T \ 

-to 

\Ri) 



(A14) 



where 



^rt f vtt\ — l ]\/[''"'' x^''' X''*(X**)~^X*^ 



(A15) 
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Prom Eqs. (A12)-(A15), the explicit expressions for M^l^j. matrices are given by 



"0 



UO+U+ 



up 

\ 1M)+U- 



UQ+U+ 



M: 



tr 



UO+U+ 



Uo+U- 



U-—UO 
Uo+U- / 



M: 



m: 



rt 



Uq—U^U- 



(uo+n+)(wo+«-) 
j uo(u+-u-) 
(uo+"+){mo+w- 

2u+ 2m_ 



(«()+«+ )(*»0+'"-) 

Uq—U+U- 

(uo+u+){uo+u-) 



UO+U+ 

A 2^+ 



«()+«- 
• 2u_ 



(A16) 



where the ij-elements of M^* matrices are defined by Fig. 10, and the subscript n — > r 
represents the incidence from a nonmagnetic medium into a resonant magnetic one. Prom 
the definition of M*'« matrices in Eq. (A14), i?i5(ki) and Tj°^(ki) correspond to M^^r and 
M^_^j., respectively. For the time-reversed waves incident with vector (— k/), scattering angle 
9f, and polarization u, M^'i^j.{—kf) matrices are same as the case of (kj,/^) but replacing i 
by (— i) in Eq. (A16), i.e., 



M^l^{-kf) = M^l^iki, z^-i), ipq = tt, tr, rt, rr), 
Tj'J{-kf) = T^'J{k,,i^ -i), <)(-k^) = i?S(k,;i ^ -i). 



(A17) 



For completeness, let us now consider the reverse case where a wave is incident "from" 
a magnetic (anisotropic) medium with Xap — Xi^ap + Xaj "into" a nonmagnetic one with 
Xa/3 — Xo^ap- The explicit forms of M^^^ matrices can be evaluated by starting with 
reversing both sides in Eq. (A12) and representing the waves as Tq — {E^J-\ E^^^), Rq — 
Ti = {Eoa,Eo^), and R, = {Er„Er^) in Eq. (A13). Then, M^l^ matrices can 
be obtained straightforwardly by 



m: 



M. 



rt 



(A18) 



where the subscript r ^ n denotes the incidence from a resonant magnetic medium into 
a nonmagnetic one. In the same way as in Eq. (A17), the M^^„(— k/) matrices for the 
time-reversed waves can be also obtained by replacing i by (— i) in Eq. (A18). 
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Finally, let us also consider the magnetic-magnetic (resonant-resonant) interface between 
upper (xa/3,up = Xup(^a/3 + Xa^up) ^nd lower {xai3,dw = XdwSap + Xaldw) resonant magnetic 
layers. By employing the 4x4 matrices involving resonant magnetic medium to both sides 
of Eq. (A12), the explicit expressions of M^^^ can be given by 



tr 



Ml 



/ 2^ 

„dw , , up 







2m5 







rt 



-m: 



tr 



(A19) 



where u'^''^^ = 

~l~ Xup,dw =•= -^up,dw -^up,dw '^^s defined in Eq. (A4). Note that these 
M*'^^ matrices for the magnetic-magnetic interfaces are apphcable to the nonmagnetic- 
nonmagnetic (nonresonant-nonresonant) interfaces simply by setting B'^^ ^^ to be zero. 



APPENDIX B: EVALUATION OF THE MATRIX ELEMENTS INVOLVING x^°^ 



To evaluate the matrix element in Eqs. (4.15) and (5.12), we assume that Xo = in Eq. 
(4.7), i.e., the first nonmagnetic medium is vacuum. Then the matrix element in Eq. (4.15) 
can be evaluated from Eqs. (4.3) and (4.10) as 



3 

X I]e*„(xi5a/3 + Xa/3)eM/3 



X l^TTJ-) Tr2^Sa(xA/3 + Xa/3)eM/3- 

j '^fzyJ) l^iz a/3 



(Bl) 



In order to evaluate the explicit expression for the above equation, let us now consider 
the case where the incidence angle 9i is small and M || x, as discussed in Appendix A. In 
this and k^fzU) — koU±, where the upper and lower signs correspond 

to j = 1 and 2, respectively, and kiz — —koUQ. Prom Eqs. (4.2), (A4), and (A7), the 
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polarization-dependent terms are evaluated by 



E e*jaiXiSap + Xal^nP = Xi(e* " 6^) + (e* X e^)xXm\\st 



a/3 



(2) 



u 



u 



0' 



for j = = a 



(B2) 



Xi - «X^^^ = li?. - ul, for j = 2, /X = (T 
-^Xi + X^^^ = -^(^*+ - ^*o)> for J = 1, = TT 

«Xl + X^^^ = - l^o)' fo^^ i = 2, = TT 

where = when xo = in Eq. (AlO). The explicit form of 2 x 2 matrix rj°^(— k/) can 
be obtained from M^_^^ in Eq. (A16) by replacing i by {—i). Then, the matrix element in 
Eq. (4.15) can be expressed by 2 x 2 matrix in terms of the polarizations of incident and 
outgoing beams, ii and as follows: 

2mo 



kl < -kJ,i/|x^°^|E;(r) > = iAkl5k^^kJk^^kf 

( 



X 



/co(m+ + mo)(m- + «o) 



= 2iAki;,R';^^{ki)Ski^kfJk,ykfy, 



(B3) 



where i?^^)(kj) corresponds to M^^^^ in Eq. (A16). Without loss of generality the final result 
in Eq. (B3) is apphcable for the case with Xo 7^ although the calculation for 2; > should 
be included in Eqs. (B1)-(B3). 

For the transmission coefficient, the matrix element in Eq. (5.12) for Xo = can be also 
evaluated from Eqs. (5.10) and (5.11) as 

E<a(Xi5a/3 + xS)e,'/3 dze-'^-'f^-''-^^'^> , 

_ o J—00 



X 



a/3 



— iAk^Ski^kf^^kiykfy 

^ rg)(-k,)7j°j(k,) ^ w , ^ (2), 



(B4) 



where the vector field E(kj,/x) in Eq. (5.10) has been used for the state |kj,/x > instead of 
the "pure" incoming wave E^(r) in Eq. (4.3). Similarly to the reflection coefficient in Eqs. 
(B1)-(B3), the matrix element in Eq. (5.12) can be expressed by a 2 x 2 matrix in terms of 
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the polarizations of incident and transmitted beams, // and j, as follows: 

A ( u "Q u 

^ _uT oIyWIIc. 1/ > — j Ak'^fi, 1 Si. I. +U++UO +U++UO 



U—+UD U-.-\-UO 

up —iUQ 



up lUp 
U-+Up U-+Up 

= iiAkiij)Tjl\]^)Sk^^kJk.,k,,, (B5) 

where Tj°'*(kj) corresponds to in Eq. (A16). Again, the final result in Eq. (B5) is 

applicable for the case with Xo 7^ without loss of generality. 



APPENDIX C: EXPLICIT EXPRESSIONS FOR ROUGH-INTERFACE Mpi MA- 
TRICES 

For the interface between upper nonmagnetic {Xai3 = Xo^ais) and lower resonant magnetic 
(Xa/3 = Xi^ai3 + Xa/j) layers, the explicit expressions of the rough-interface M^'i^^ matrices 
can be given by 



fl- V. 



(CI) 



where, from Eqs. (5.6) and (5.19), 



i-v. 



I -v. 



(m++«o)^ 



+ 



d\+^+d, 

(u++«o)^ 



(+) 



£,( + )_ 



D 



{+) 



(«_+uo)^ 



(-) 



+ 



D 



(-) 



(-) 



D 



(-) 



(-)^^(-) 



(-) 



D 



{-) 

4 

!2 



(+)^n(+) 



(•it++-U.())2 



(+) 



d\+^+d._ 

(U++M0)^ 



+ 



D 



-) 

(«-+«o)^ / 



(+)_ 



2 

.-"■0 
~~2 



\ 



and 



n— 



— Z 



. "(1 g "(1 . 



(C2) 



(±) 



D 



(±) 



(Xi 
(Xi 



- Xo)e" 
-Xo)e 



D 



(±) 



(C3) 
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Here, {xi — Xo) — {^X + uV)/'^ — ul and B' — {u^ — u'i_)/2 can be used from uq — 6f + Xg 
and u± = ^Jef + Xi ± B' . 

For the reversed interface between upper magnetic (resonant) and lower nonmagnetic 
layers, -M^^^ matrices can be also given by 



Mi: 



(I-V,^„)-^(R(°)„ + U,^„), 



(C4) 



where 



I- V. 



(u++wo)^ 

r,(+)_n{+) I ' 

I-V I-V =l-V 



■yCO) — /I/** — 



2U4 



2it_ 



2u 



±_ 



• 2u_ 



(C5) 



\ tt++llo 



In the same way as Eq. (A18), other M^'^ matrices can be given by 



(C6) 



Finally, for the magnetic-magnetic (resonant-resonant) interface between upper resonant 
magnetic {x7(3 = x7^ap + Xai""^) and lower resonant magnetic (x^^ = xf^^ap + Xaj''^") 
layers, Mf^^ matrices can be also given by 



Ml 



where 



I-V, 



(I-V,^,)-^(R(°l, + U,^,), 
(I-VU,)-^T(°)„ 



(C7) 







(udw)2_(„;jP)2 





±11 £is 





D\-'^-Di 



(«i^)2-(«:ip)2 y 



(C8) 



and 



(±) 



(±) 



(xf - x^^)e-^K"±<^)^'^?, = (B^^ - B;je-^K-±<T<^i 

(Xf - x?'')e-*("^-^^"-^)^'^?, i^?) = (S^. - )e-*(«-±«-)^<^i (C9) 
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and (I — \'r-^r) corresponds to (I — Vr_»r) when switching the upper and lower layers, and 
Ti% corresponds to M^^^ in Eq. (A19). Here, (xi^-xT) = [(u^y + (u^^y]/2-[(u^f + 
(M:iP)2]/2 and {B'^^ - B'^^) = [(wf )2 - (m^)2]/2 - [{u'^f - {uIf]/2 can be used from 
^up.dw _ _|_ ■y^''^"" J- i?^p^j^. In the same way as Eq. (C6), two other matrices M*!^^ and 
-^r^r t)e also obtained from M^^^ and M*^^^ in Eq. (C7), respectively, by switching the 
upper and lower layers. We should mention again that these rough-interface M^'^ matrices 
for the magnetic-magnetic (resonant-resonant) interfaces can be reduced to the cases for 
the nonmagnetic-nonmagnetic (nonresonant-nonresonant) interfaces by setting B'^^ ^^ to be 
zero. 



APPENDIX D: SOLUTIONS OF SELF-CONSISTENT MATRIX EQUATIONS 
FOR NONMAGNETIC INTERFACES 

For nonmagnetic interfaces (|M| = 0) and cr — > cr polarization, simply 

u^^u_^^Jej + xi = 141 Ao. (Di) 

Inserting this in Eq. (A16) modified for (— kj) and using Xi ~ Xo = (I^P ~ l^2p)/^o> 
self-consistent solution for the reflection coefficient (k^ = k[ and Oi — Of) in Eq. (5.6) can 
be reduced to a scalar as 

R={1-V)-\R^^^ + U), (D2) 

where 

U = l^^l^l^^l (e-U\>^^\+\>''^\?-l - l) , (D3) 
and fcj^ = I A;* I, and kiz = —k^^ = —\kz\- Then, we obtain 

\kz\ + l^zl 

which is consistent with the Nevot-Croce form.^^ 







\K\ 


\kz\ 
\kz\ 


+ 


\K\ 
\K\ 


1 kz 1 ~l~ 1 k^z 1 
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Similarly, the self-consistent solution for the transmission coefficient in Eq. (5.19) can be 
reduced into a scalar as 



where 



Then, 



i=l,2 j=l,2 



y> = l_ f^-\{\kz\-\kl\fal^ 

T'S = rr = -Jy^. (D6) 



2\k 

which is consistent with the Vidal- Vincent form.^^ 



\K\ + \ki\ ' ^ ' 



APPENDIX E: RECURSIVE 2x2 MATRIX FORMULAE FOR MULTIPLE IN- 
TERFACES 

For multiple interfaces, additional phase differences between different interfaces should 
be taken into account to extend the results for a single interface in Appendix A. Following 
Ref. 15, MnXi matrices for the n-th interface between n- and (n-l-l)-th layers can be modified 
from Eq. (A15) as 

M*Vi = M«F-i, M*;l = M*^ M:,\, = F-'M^^-'F-\ M:1, = F-'M^\ (El) 

where M^* are the 2x2 matrices obtained for a single smooth interface in Appendix A, 
depending on whether the upper and lower layers on the n-th interface are nonmagnetic or 

magnetic ones, respectively, 

^-ikou+,„d„ g 



Fn 



(E2) 



and u±^n and dn represent the refracted angle defined in Eq. (A7) and the thickness of the 
n-th (upper) layer, respectively. For nonmagnetic layers, u±^n reduces to iio,n in Eq. (AlO). 
Rn and Tn are the vectors {Rn,i,Rn,2) and (T„_i,T'„^2) representing the two waves reflected 
and transmitted, respectively, at the top of the n-th layer. (In Ref. 15, they are defined at 
the bottom of the n-th layer.) 
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Introducing VF^* matrices following Ref. 15, which are defined by 

To 




n ' ' n 



(E3) 



and using the recursion formulae involving Mf^^ matrices at the n-th interface, i.e.. 




T 



(E4) 



■n+l 



yields the following recursion formulae for FF^^ matrices: 



rt 

n+l 



(E5) 



where and Bn are defined by 



5n = Vr7(l-M;*+,Pr*'-)"\ 



(E6) 



Here WJ^ determines the reflectivity of the whole multilayer, Rq = W^Tq {Rn = 0), from 
Eq. (E3). 

Finally, the field amplitudes T„, R„ inside the layers can be obtained from Eqs. (E3)-(E6) 

by 



R„ = {l-M:x^w^^y\M:l,Rr.+^+M:x,wyo), 



(E7) 



which must be progressively applied to all the layers starting at the multilayer substrate 
where Rn — 0. 
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FIG. 1: Schematic of scattering geometry and sketch of the chemical (or structural) {zc{x, y)) and 
magnetic {zm{x,y)) interfaces, which can be separated from one another by an average amount A. 
Grazing angles of incidence {0i) and scattering {Of), the wave vectors kj and k/, and the photon 
polarization vectors of incidence (e^=o-,7r) and scattering (e,^=cr,r) are illustrated. Small arrows 
represent the possible orientations of the magnetic moments around magnetic interfaces. 




FIG. 2: Schematic of an ideal interface with undisturbed states E(ki) and E (— k^). Note two 
possible waves for each of the reflected and transmitted wave vectors. 
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FIG. 3: Calculated x-ray resonant magnetic reflectivities at the Gd Ls-edge (7243 cV) from Gd 
surfaces with different interfacial widths for structural (ctc) and magnetic ((7^) interfaces: (a)-(c) 
(Jc = SA, am = 3A. (d)-(f) (Tc = 3A, am = SA. (g)-(i) same as (d)-(f), but with a 20A magnetically 
dead layer. Top panel: reflected intensities of the a ^ a (solid lines) and a ^ tt (dashed lines) 
channels, and the differences between the reflected intensities for right- (/+) and left- (/_) circularly 
polarized incident beams (circles). Middle panel: Natural logarithms of the reflectivities with 
interface roughnesses normalized to those from ideal systems without roughness as a function of 
the square of the wave-vector transfer. Solid, dashed, and dot-dashed lines represent a ^ a and 
a ^ TT scattering, and the differences between /+ and respectively. Bottom panel: Normalized 
scattering density profiles for charge (solid lines) and magnetic (dashed lines) scattering. 
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FIG. 4: Charge (a) and magnetic (b) x-ray scattering amplitudes, fc^m around the Gd L2-edge 
obtained from the absorption measurements for a [Gd(51 A)/Fe(34 A)]i5 multilayer. The vertical 
lines indicate the photon energies, where the x-ray resonant magnetic reflectivities in Fig. 5 were 
calculated. 
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FIG. 5: Calculated x-ray resonant magnetic reflectivities from a [Gd(5lA)/Fe(34A)]i5 multilayer 
for different incident photon energies indicated in Fig. 4: (a) 7926 eV, (b) 7929 eV, (c) 7931 eV, 
and (d) 7935 eV. Both structural (charge) and magnetic interface roughnesses are ac,m = 4.7 A and 
3.6 A for Fe/Gd and Gd/Fe interfaces, respectively. The solid lines represent (/+ + /_) intensities 
and open (filled) circles represent the positive (negative) values of (/+ — I-) intensities. 
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FIG. 6: Calculations with different magnetic interface roughnesses: (a) and (c) am = 2.1 A , and 
(b) and (d) am = 6.2A. All other parameters and symbols are same as those in Fig. 5. 
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FIG. 7: Models of magnetic structures in Gd layers. Uniform magnetization (A), ferromagnetic 
moments only near the Gd/Fc interfaces (B), and ferromagnetic moments near the centers of 
Gd layers between magnetically dead layers (C). While interfaces with "cTc, (7^" represent both 
structurally and magnetically mixed interfaces, interfaces with "ctc" (or "cr^") represent purely 
structural (or magnetic) interfaces. 
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FIG. 8: Calculated (/+ — /_) intensities for different magnetization depth profiles in Gd layers. 
In (a)-(c) ferromagnetic layers exist only near the Gd/Fe interfaces [Fig. 7(B)], and their layer 
thicknesses are 4.6 A (a), 8.4 A (b), and 12.8 A (c). In (d)-(f) ferromagnetic layers exist in the 
middle of Gd layers and are sandwiched between magnetically dead layers [Fig. 7(C)], and the 
layer thicknesses of the dead layers are 4.6 A (d), 8.4 A (e), and 12.8 A (f). All magnetic roughness 
amplitudes are am = 4.2 A, which is effectively same as cJc, and the photon energy is E = 7929 
eV. All other parameters and symbols are same as those in Fig. 5. 
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FIG. 9: (/+ + /_) [(a)] and (/+-/_) [(b)] intensities measured (symbols) from a Fe(34 A)/[Gd(51 
A)/Fe(34 A)]i5 multilayer near the Gd L2-edge (7929 eV). The lines represent the best theoretical 
fits with the model (B) in Fig. 7. Note that the colors of symbols and lines in (/_(_ — /-) intensities 
are different for opposite signs of the intensities. 
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FIG. 10: The representation chosen for the elements of M^'^ matrices with the polarization bases 
of the incident and reflected (or transmitted) waves. The polarization basis is given by (eo-, e^r), as 
shown in Fig. 1, for the waves in the nonmagnetic medium and (e^^^ e^^^), as defined in Appendix 

A, for those in the resonant magnetic medium, respectively. 
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